Evolution of domain wall networks: the PRS algorithm 
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The Press-Ryden-Spergel (PRS) algorithm is a modification to the field theory equations of mo- 
tion, parametrized by two parameters (a and implemented in numerical simulations of cosmo- 
logical domain wall networks, in order to ensure a fixed comoving resolution. In this paper we 
explicitly demonstrate that the PRS algorithm provides the correct domain wall dynamics in A'^ + 1- 
dimensional Friedmann-Robertson- Walker (FRW) universes li a -\- jS /2 = N , fully validating its use 
in numerical studies of cosmic domain evolution. We further show that this result is valid for generic 
thin featureless domain walls, independently of the Lagrangian of the model. 
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I. INTRODUCTION 



The dynamics of cosmological domain walls has been 
investigated using both high-resolution numerical simula- 
tions and a semi-analitical velocity-dependent one-scale 
(VOS) model [l|-l7|. Most of these studies were moti- 
vated by the suggestion Q that a frozen domain wall 
network could be responsible for the observed accelera- 
tion of the Universe (see also P-O^ ) . Although, current 
observational constraints on the equation of state param- 
eter of dark energy strongly disfavor domain walls as a 
single dark energy component [l^ [ij] , they are unable to 
rule out a substantial impact of a frustrated domain wall 
network on the acceleration of the Universe around the 
present time. However, analytical and numerical results 
strongly support the conjecture that no frustrated do- 
main wall network, accounting for a significant fraction 
of the energy density of the Universe today, could have 
emerged from realistic phase transitions. These results, 
on their own, seem to rule out any significant contribu- 
tion of domain walls to the dark energy budget. How- 
ever, they rely heavily on the validity of the so-called 
Press-Ryden-Spergel (PRS) algorithm used in cosmolog- 
ical domain wall network simulations. 

Domain walls have a constant physical thickness and, 
consequently, their comoving thickness decreases propor- 
tionally to the inverse of the cosmological scale factor. 
In numerical studies of cosmological domain wall evo- 
lution the rapid decrease of the comoving domain wall 
thickness would be serious problem since it would im- 
ply that domain walls could only be resolved during a 
small fraction of the simulation dynamical range. The 
PRS algorithm is a modification to the field theory equa- 
tions of motion, implemented in numerical simulations of 
cosmological domain wall evolution, allowing for a fixed 
comoving resolution. It has been argued that the PRS 
algorithm [l[ , provides the correct domain wall dynamics 
in 3 -f 1 dimensions, as long as a + (3 /2 = Z {a and /3 



are the PRS algorithm parameters of ref. jj ) . Although 
this claim is strongly supported by numerical tests it has 
never been proven that the same Nambu-Goto effective 
action is recovered in the thin wall limit. In this paper 
we eliminate this shortcoming, extending the analysis to 
generic thin featureless domain walls in FRW universes 
with an arbitrary number of spatial dimensions. 



II. THE PRS ALGORITHM I 

Consider the Goldstone model with a single real scalar 
field (j) described by the Lagrangian 



c = x- v{4>) 



(1) 



where X — ~(f>^^(f)'^^ /2 and V{(f) is the potential. This 
model admits domain wall solutions if the potential, 
V{4>), has, at least, two discrete degenerate minima. 
Varying the action. 



S = 

with respect to the scalar field, 
ing equation of motion 



(2) 



one obtains the foUow- 



dV 



(3) 



Here g — det(f;ct^) and gap is the metric tensor. In this 
paper the Einstein summation convention will only be 
used with greek indices (such as in Eq. ([3])). 
In a flat FRW universe, the line element is 



ds^ = a'^{vi){-dTf + dx • dx) , 



(4) 



where a(r/) is the scale factor, 77 = dt/a is the conformal 
time, t is the physical time and x are comoving coordi- 
nates. The equation of motion for the scalar field 4> given 
by Eq. ([3]) becomes 



,dV 



(5) 
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where a dot represents a derivative with respect to con- 
formal time, TV is the number of spatial dimensions, 
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H = a/a and is the comoving Laplacian. A static 
straight domain wall solution oriented along the x di- 
rection can be obtained by choosing initial conditions 
such that ip = (l){x) with ^ = and cj) — (we take 
X = (xi, X2,X3) and xi = x). Eq. ([5]) preserves the phys- 
ical thickness of the domain walls so that the comoving 
thickness is proportional to a^^. This is a problem for 
cosmological domain wall network simulations since the 
comoving thickness of the domain walls decreases very 
rapidly and can only be resolved during a small fraction 
of the simulation dynamical range. 

The PRS algorithm consists of the following modifica- 
tion to the equations of motion 



dV 
d(j) 



(6) 



where a and (3 are constants. By taking /3 = it is pos- 
sible to fix the comoving thickness of the domain walls so 
that they can be resolved throughout the full dynamical 
range of the simulations. Moreover, it was shown that 
ifa + /3/2 = 3 the dynamics of planar domain wall in 
a 3 -|- 1-dimcnsional FRW universe would be maintained 



a- 



a rapid change of occurs only in the direction orthog- 
onal to the wall [l^- It is convenient to choose spatial 
coordinates (u, w, z) such that locally the walls are coor- 
dinate surfaces satisfying the condition u — constant. In 
this case, the domain wall is parameterized by the coor- 
dinates w and z and it moves along the it-direction. It 
is useful to choose an orthogonal coordinate system in 
which w = constant and z = constant are lines of cur- 
vature so that the coordinate curves coincide with the 
principal directions of curvature of the surface defined 
by u = constant. It is always possible to construct such 
a coordinate system in the vicinity of any non-umbilic 
point (in which the two principal curvatures exist and 
are not equal) of a surface embedded in a flat space [11] . 

If the domain wall has velocity v, then the domain wall 
solution is still be given by </> = 4>s{l) with 



dl 



— -7 -^--^-0 (11) 



III. THE PRS ALGORITHM II 



where dsi — \dfi\ is the arc length along direction Ui and 
dfi — hiduiUi (ui is the unit vector along the direction 
Ui). We shall use the gauge freedom to choose a coordi- 
nate u which measures the arc-length along the direction 
perpendicular to the domain wall, so that = 1 and 
dsu = du. 

Therefore, one has 



Changing the space-time coordinates, in Eq. ([6]), from 
(7?,x) to (^,y), defined by 





d 


1 d 

al^/^dT]' 


(7) 




y 


= a'^/^x, 


(8) 


de ^ 






(9) 



yields 



where = a'^^l and H = a-^/^H. 

In Minkowski space-time (a = 1) a planar static do- 
main wall solution oriented along the y direction will be 
given hy (j) = (j}s{l) with 



dV 



dP 



(10) 



with I = y {we take y = (2/1,2/2,2/3) and yi = y). If the 
domain wall is boosted along the positive y direction, the 
planar domain wall solution to eq. Q is still (j) = 0s (0 
but now I = 7(2/ — w^) where v is the domain wall velocity 
and £, = t. In this case dl/d^ = —-fv and 01/ dy = 7. 

Consider the more general case of a curved domain wall 
in a 3-1-1 dimensional flat FRW universe. The generaliza- 
tion to A'^-l-l dimensions is trivial and for simplicity it will 
be made only at the end of the section. In the following 
we shall assume that the thickness of the domain walls 
is very small compared to their curvature radii, so that 



dl 
d\ 



dP 



dcf) d4>s 

_ d{jv) d4>s 
dl 



On the other hand, taking into account that 
and /i„ = 1, 



V^0 = 



1 



1 dhw 
/i,„ du 



d_ 
du 

hy du 



hit} hz 



hu du 
90 9^ 
du du^ 



(12) 
(13) 

' = 0(C,'") 

= (14) 



The curvature of a curve parameterized by p is defined 



as — 



kpl where kp = dep/dsp, Gp is the unitary tan- 



gent vector to the curve and dsp is the arc-length. The 
principal curvatures of a surface, defined by a constant 
u = uq, are given by — k^, • u and A:^ = k^ • u with 



k?i) — 
k. = 



1 /9w 

i_ rdz- 

hz \dz 



(15) 
(16) 



The vectors u, w and z form an orthonormal but, in 
general, non-coordinate basis. Their derivatives can be 
calculated using 



duj 
duj 



1 dhj . 
hi dui 



1 dhi 
^ hk duk 



Ufc • 



(17) 
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Hence, 



1 dhu 
du 



■u , 



1 - ^1Q^ 



The relevant curvature for domain wall dynamics is the 
extrinsic curvature, i.e. the "bending" of the wall in rela- 
tion to the flat embedding universe. Mathematically this 
is measured by the curvature parameter 



1 dh^, 



hw du 

Therefore, Eq. (1141) can be written as 



1 dh^ 
h, du 



du^ 



(19) 



(20) 



Inserting this in Eq. ([5]), taking into account Eqs. ((T^) 
and (fT3|) and the fact that d'^4>/du'^ — ^'^dP4>s/dl'^, one 
obtains 



dl 



dV_ 

'del) ' 



(21) 



where 



a - 



^]Ujv + JCj. (22) 



Taking into account Eq. ([TU]). we conclude that ~ 0. 

Notice, however, that dsw and dsz are not the comov- 
ing arc-lengths since the comoving space coordinates have 
been scaled by a factor a~^^^. The comoving curvature 
parameter is instead 



(23) 



Changing into the original variables (77, x), we finally find 
that 



= 0. 



(24) 



By setting the parameters a and /3 to their original values 
(a = /? = 2) we find that, if the modified equations are 
to yield the correct dynamics in a 3-|-l-dimensional FRW 
universe, we must have that a + (3/2 — 3. 

In a A^-|-l-dimensional FRW universe domain walls are 
defects with — 1 spatial dimensions whose dynamics is 
still given by Eq. Q (see ref. [l3| for an analytical study) 
with 



N-l 

1=1 



(25) 



IV. GENERIC DOMAIN WALL MODELS 

In this section we show that the main result of the pre- 
vious section (Eq. describes the dynamics of generic 
thin domain walls, independently of the Lagrangian, 
C{(j),X), of the model. We will follow closely the deriva- 
tion presented in ref. [6] where the validity of Eq. ([M)) 
has been demonstrated for planar domain walls. Varying 
the action, 



S 



dt 



d^Xy^C{(l),X) 



with respect to (f), one obtains 



'-9 



-c., 



(26) 



(27) 



where C^x — dC/dX and — dC/d<j). 

Assuming a -I- 1-dimensional FRW metric given by 
Eq. (HI and the transformations given by Eqs. ([7]) and 
dHl) one obtains 



d_( 



Vy0-£A-V2 = /:,^(28) 



with a = N — 1 and (3 = 2. Notice that, in Minkowski 
spacetime, a planar static domain wall solution oriented 
along the y direction will be given by (j) = (j^sil) with 



^ — ( r ^ 
dl \ dl 



(29) 



with / = y. 

Consider the coordinate system {u, w, z) as described 
in section HIl Suppose that the wall is moving along the 
direction u with velocity v. Taking into account that 



Vv/: 



dC 



du du 



(30) 



as well as Eqs. (|lim3p and (PP)) . the equation of motion 
yields 



dl \ dl 



dl 



(31) 



Again, since (f){l) must be a solution of Eq. (P^ one has 
J- = Q and, consequently, Eq. (j24p remains valid. Fur- 
thermore, although only models with a single real scalar 
field have been considered, it is straightforward to ver- 
ify that Eq. ([M]) describes the correct thin domain wall 
dynamics in the context of generic models with various 
scalar fields. 



Here, are the curvature vectors associated with the 
N — 1 coordinate curves of the domain wall. Hence, the 
dynamics of thin domain walls is unaffected by the PRS 
algorithm as long as a+(3/2 = N (the original parameters 
were a = — 1 and /3 = 2). In particular, the dynamics 



planar (k = 0) domain walls is such that W7 cx a 

„-N 



-Q-/3/2 



V. DOMAIN WALL DYNAMICS IN 2 + 1 
DIMENSIONS 

The world history of an infinitely thin domain wall 
in a flat FRW universe can be represented by a two- 
dimensional world-sheet with x = x(?7, cr), obeying the 
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usual Goto-Nambu action. The equations of motion take 
the form 

x + 2H(l-x2)x = e-i(e-ix')' (32) 
e = -2Hek^ , (33) 



with 



i-x' = 0. 



J2 \ 2 



1 -X2 



(34) 
(35) 



where dots and primes are derivatives with respect to rj 
and (T, respectively. 

Let us define unit normal and tangent vectors as 



X X 

u = - , w = — , 

V C 



(36) 



where v{ri,a) — |x| and C{ri,a) = |x'|. Eq. ((35l) can now 
be written as e = 7C with 7 = (1 — v^)^^^^. Therefore, 
the left hand side of Eq. ([5^ is given by 

k + 2-H (1 - x^) X = i-u + ufl + 2n (1 - v^) vu , (37) 

where u is proportional to w. Moreover, the right hand 
side of Eq. ((32|) gives 



7 (9s \ 7 
7^ as 



(38) 



where we have taken into account that dw/ds = ku and 
the fact that the physical length along a 2-dimensional 
domain wall is given by ds — \dx\ — Cda. Henceforth, 
the normal component of eq. ([5^ yields: 



(39) 



v + {l-v^){2nv- k) = 0, 



which confirms Eq. (|24p in the particular case with N 
2. 



VI. CONCLUSIONS 

In this paper we explicitly demonstrated that the PRS 
algorithm provides the correct dynamics of thin feature- 
less domain walls in FRW universes with an arbitrary 
number, N, of spatial dimensions, if a + /3/2 = iV. Our 
results fully justify the use of the PRS algorithm in nu- 
merical studies of cosmological domain wall network evo- 
lution. Although, fixing the comoving thickness of the 
domain walls, using the PRS algorithm, increases artifi- 
cially the impact of the junctions on the overall network 
dynamics during the course of the simulations, this effect 
is negligible for the light junctions usually considered in 
such simulations. 
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